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Abstract: 

^ ■ Evolution equations for multiplicities in QCD cascades can, both in the parton and dipole 

picture, be used to estimate corrections beyond the formal accuracy of the modified leading 
log approximation (MLLA). The differences between the two pictures, and other uncertain- 
ties beyond first order MLLA corrections, are here investigated in some detail. For example, 
I discuss how some colour suppressed corrections, which cannot be determined without bet- 
ter knowledge of non-perturbative QCD, are related to the question of colour reconnection. 
A generalized evolution equation for the dipole cascade is also presented. 
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1 Introduction 



The modified leading log approximation (MLLA) |]I], of QCD cascades systematically 
includes corrections suppressed by a relative factor ^/a^, where a^ is the coupling of the 
strong interaction. Confronting MLLA predictions with data on multiplicities in quark and 
gluon jets, iVq and A^g, higher order corrections are found to be important 0, |]. These 
corrections, in particular energy conservation effects, have been estimated by expanding the 
MLLA evolution equations beyond their order of formal reliability p, |^, |^ . Thus, no exact 
and systematic inclusion of 0{as) effects is available. 

One contributor to 0{as) corrections is the class of Feynman diagrams where the gluons 
are not strongly ordered, i.e. where two gluons have similar transverse momentum, or where 
the hardest gluon has a transverse momentum close to the kinematical limit. The iterative 
cascade formalism is not designed to include such diagrams in a proper way, and thus "mod- 
erately ordered" gluons introduce uncertainties to the 0{as) corrections, which may have 
numerical consequences. 

The iterative cascade is a branching process, where each emitter can split into two new ones, 
which continue to emit independently. However, colour suppressed terms, which reflect 
possible interferences between colour charges due to the finite number of colours, do not fit 
naturally into this iterative picture. According to [0, approximately 10% of the important 
energy conservation corrections to Ng/Nf^ are related to these uncertain effects, which cannot 
be fully determined within the perturbative QCD cascade formalism. In this paper I discuss 
their connections to the question of colour reconnection p], ^. 

A simple way to investigate the uncertainties in the estimates of higher order corrections is 
to compare two alternative evolution equations, based on the parton cascade and dipole 
cascade [|10|, respectively. These are both equivalent within the MLLA accuracy, and their 



differences are a guide to the uncertainties at higher orders. In this paper, I compare the 
theoretical results in the dipole and parton pictures more directly than has been done before, 
by expanding the dipole evolution equations in ^/Q^, in the same way as has been done for 
the parton cascade. 

I also present new dipole evolution equations, alternative to the ones in |T^. The new 
derivation is less transparent and the result is mathematically more involved but, in spite of 
such shortcomings, it is presented for two reasons: First, discussing the differences between 
the two alternative dipole evolution equations is a powerful way to examine the uncertainties 
beyond first order MLLA corrections. Second, the recoil corrections, presented in [^ for 
average multiplicities only, are here treated in a more general way, suitable for calculations 
of all multiplicity moments. 

The outline of this paper is as follows: In section 0, I present the investigated evolution 
equations. In section ^, their predictions on some multiplicity quantities are compared. The 
differences and other uncertainties at higher order corrections are discussed in section |[ In 
section ^ a relation to colour reconnection is discussed. The paper ends with a summary in 
section O. 



2 Evolution Equations 

In this section, I list three alternative evolution equations for the multiplicity distribution in 
QCD cascades. Though this paper focuses on average multiplicities, I will when possible give 
the evolution equations for the more general generating functions. The generating function 
to a distribution, P(A), of some quantity A, is defined as 

V{z)= j dXP{X)(l + zY. (1) 

While the distribution for two indepently contributing sources is described by a convolution, 
the corresponding generating function is simply multiplicative, 

P,+,(A) = J dX^dX2Pa{Xl)Pb{X2)S{X - Ai - A2) => Va+b{z) = Va{z)V,{z). (2) 

(In case of a discrete variable, e.g. multiplicity, the integrals are replaced by sums, and the 5 
distribution by a Kronecker 5.) Expanding V in powers of z gives different moments of the 
multiphcity distribution, and thus knowledge of V gives full knowledge of the distributions. 



2.1 Part on Equation 

The evolution equations for generating functions for the multiplicity distributions are in the 
parton cascade picture |^ given by 



V'g{1) = JdxKS{x)2a[VG{l + Hx))VG{l + Hl-x))-rG{l)] 

■ J dxK^{x)2a[VF{l + \n{x))VF{l + ln(l - x)) - Vg{1)], (3a) 



2N, 



V'pil) = jf f dxK^{x)2a[VG{l + Hx))Vf{1 + ln(l - x)) - Vf{1)]. (3b) 

Here Vf represents a quark jet and Vg a gluon jet. x is the momentum fraction of the emitted 
gluon (or quark, in case of g— >qq splitting), Uf is the number of active quark flavours, N^ = 3 
is the number of colours and Cp = {Nc/2){1 — N~'^) = 4/3. The dependence on the auxiliary 
variable z is suppressed in the notation, and as kinematical variable is used / = ln(Q/A), 
where Q is the virtuality of the jet and A is the scale parameter of the perturbative cascade. 
The running coupling, depending on the transverse momentum k±, is written as 

^ = ^ = ?^(1_?AH^)^0(.-). (4) 



where 



2n j3qk PqK 



K = \n{kl/A^), (5) 

(3o = l (llA^c - 2nf) , (3i = ^ {UN^ - n^ibN, + SC^)) . (6) 



The splitting kernels K can be written 



X X 

Kq must be chosen so that the sum 

\[K'i{x) + K^{1 - x)] = i^ + ^ + x(l - x) (8) 

agrees with the Altarelli-Parisi splitting function for g^gg. The choice made here is conve- 
nient for comparison with the dipole formalism. 



2.2 Dipole Equation 



In [g, the dipole evolution equations presented in [^ are supplemented with recoil correc- 
tions. Only the average effect of recoils is considered, and it is therefore preferable to use 
the equations for the average multiplicity only, and not higher moments. The result is 

A^g (L) = a{L - Cg) [l - cfe)a(L - Cg)] iVg(L - Cg) 

K(L) = ^a(L-Cq)[l-4^)a(L-Cq)]iVg(L-Cq). (9b) 

Here Ng{Nq) is the mean multiplicity of a gluon(quark) jet, defined as one hemisphere of a gg 
(qq) system stemming from a point source, with invariant mass y/s. This scale is represented 
by the logarithmic variable 

L = \n{s/A^). (10) 

The number Cg(cq) represents corrections related to the non-singular part of the gluon emis- 
sion kernel Kq{Kp), in the following referred to as "kernel corrections". The number 
c[^\cl.^^) is related to recoil effects in a gg(qq) dipole, while Cp reflects energy conservation 
in pair production g— *qq. The values are 

c = 11 c = ^ c = 1^ (11) 

cfe)=2(--^) c(^) = -AX_^ n2) 



2.3 Generalized Dipole Equation 

In appendix, I derive dipole evolution equations alternative to Eqs (|^). The recoil corrections 
are here treated in way suitable for calculations of all multiplicity moments. I therefore refer 
to these formulas as the generalized dipole evolution equations. The equations are 



[ln(Pg(L + Cg))]' = i5(L) + a(L)Pg(L + Cg)exp(-Cgi3(L))^i?fe)i3"(L) + 



[ln(P,(L + c,))]' 
B'{L) 



+ 



n{a{L) 



n=\ 



^-(^+"^)exp[(c,-^c,)i3(L)]E/?iP)S"(L)-l 






A^. 



n=0 



B{L) + a(L)Pg(L + cg) exp(-CgS(L)) ^ i?i^)S"(L) 



n=l 



a(L)[Pg(L + Cg)exp(-Cgi3(L))-l]. 



(13a) 

(13b) 
(13c) 



Here B is an auxiliary generating function, whose moments in the generating function variable 
z can be eliminated, replaced by moments of the physically more relevant functions Pg and 
Pq, which represent the multiplicity in a gg and qq dipole, respectively. Note that a gg 
system consists of two gg dipoles, one from each triplet charge to the matching anti-triplet 
charge of the other gluon. The numbers R^^ ■, defined in Eqs (|5P|), (^) and (^), reflect 
energy conservation. In particular, we have c^f' = —Ri. 



3 Comparison 



Both the parton and dipole pictures have been used to predict multiplicity observables. In [Q 
the generating functions of Eqs (H) are expanded V{l + e) = V{1) +eV'{l) + ..., to get parton 
cascade evolution equations expanded in y/^. In 0, the dipole evolution equations, Eqs (H), 
are combined to 



N'JL + Cg - Cq 



Cf 



.(q) 



[l-{4^^-4-^)a{L)]N'(L). 



(14) 



Using the experimentally well determined A^q as input, and a boundary condition scale Lq 
where N^^Lq) = Ng{Lo), this equation is used to predict A^g for scales L > Lq. These parton 
and dipole approaches are confronted with data in 0]. 

In this section, all evolution equations are expanded in y^ in the same way, thus enabling 
a more direct comparison of the cascade pictures. To explore the differences, it is enough to 
consider average multiplicities, with corrections up to second order in \/^. 



With numbers Oj and Vi defined as in [^ from 






y/E(l- 2aiy/E - 4a2a + C(a^/^ 
ro (l - 2riv^ - 4r2a + 0{a^^^ 



(15) 
(16) 



the evolution equations for average multiplicities can be written in the following form, cov- 
ering corrections up to second order. 





parton 


dipole 


generalized dipole 


%2 

4g2 


^V2 

-4f8 - 4 — 

^0 


1-^2 _ „{g) 
2^q ^r 


2 IS S'' ^r *^3A^3 



Table 1: Table with second order coefficients which differ between the considered cascades. 
The values of the quantities are given in the text. Their interpretation is discussed in sec- 
tion 1^. Second order terms common to all cascades are given in Eq f|7?|j. 



rif 



3 

4 
5 



0.2 

parton dipole gen. dip. 



-0.379 -0.240 -0.015 
-0.339 -0.200 0.030 
-0.300 -0.161 0.074 



^2 

parton dipole gen. dip. 



0.426 0.620 0.633 
0.468 0.663 0.680 
0.510 0.705 0.727 



Table 2: Numerical differences on second order corrections to the anomalous dimension and 
the multiplicity ratio A^g/A^q. 



n: 



N, 



Cr' "" 



F 



^N,-ic, + ^){aNj + A 
aNg - CqiaNgY + Aq2a'^Ng. 



roSN, 



+ 92 



a^N, 



(17a) 
(17b) 



The different cascade expressions for g2 and g2 are shown in Table |l|. The coefficients 
defined in [^, reflect recoil effects in the parton cascade. The values are 



V2 



67_7r2 
36 ~Y 
Matching orders of a/^ in Eqs (|T7|), we get 



«8 = -- 



7 



«10 



71 



ro 



Cf' '" 2 ^6iV3 8iV,^ 



ri 



Cg Cq vif 



2a2 

r2 






Tii 



8Nr 



ANr 



ro3iV, 



in 



QNf^ 



92, 



Cq nf 

'^2 '^ ro3N, 



Ji -—)-92 + q2- 



(19) 

(20) 
(21) 



The numerical impact of the equation differences to the anomalous dimension N' /Ng and 
multiplicity ratio Ng/N^ is presented in Table |^ and Fig |l|. We conclude that there is a 
considerable difference between the parton and dipole pictures, also when precisely the same 
procedure has been invoked to reach predictions on physical observables. The uncertainties 
represented by differences in 02 and r2 are larger than the third order corrections calculated 
for the parton formalism ||^. 
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Figure 1: Multiplicity ratios to different orders in yci- The chosen example is rif = A and 
A = 0.22GeV. d is given by Eq (^. left: Anomalous dimension N'^/Ng as \f^ (solid line), 
■\/^(l — 2aiv^) (crossed) and y/^{l — 2aiy/^—4a2<y) for the parton equation (dashed), dipole 
equation (dotted) and generalized dipole equation (dash- dotted), right: Ng/N^ as Tq (solid 
line), ro(l — 2riv^) (crossed), and rQ{l — 2ri\^ — 4:r2a) for the parton picture (dashed) and 
dipole picture (dotted). 



4 Discussion of Corrections up to Second Order 

In this section, I discuss the different terms in Eqs (|17|) and Table [I|, and also some not 
included anywhere yet. Many corrections can be visualized with the logarithmic phase space 
triangle. Following the colour ffow connecting the partons, we label the original ones 1 and 
3, and the emitted gluon 2, and let Xi be the scaled energy 2Eij ^ of parton i. With 
kinematical variables 



V 



In 



X\ 



2 VI -X3 
the constraint 2A;_l cosh(?/) < ^i implies 



fc 



s(l -a;i)(l-X3) 



\V\< 1/max 






\-4k\ls^ 



-XL 



K 



(22) 



(23) 



Thus the allowed phase space for gluon emission is approximately a triangular region in 
the y, K-plane, cf. Fig ^a. As an introduction to this phase space picture, I discuss how it 
illustrates leading order results and first order corrections, thereafter continuing with second 
order corrections of increasing uncertainty. In the following, the figure represents a gg dipole. 



4.1 Simpler Terms 



After the emission of a gluon at yi,Ki, the additional phase space for further emissions can 
be represented by a fold of height ki, as in Fig ^. The fold contributes with a multiplicity 



K=ln(kl/A^) 




Figure 2: a) The phase space for gluon is, to leading order, a triangluar region in the y, k 
plane. A gluon emitted atyi, ni enlarges the phase space with a double-sided fold of height ki. 
This implies that the triangular phase space represents the leading order relation N'' = aNg. 
b) Illustration of some subleading corrections: The dark strips illustrate effective phase space 
reductions representing kernel corrections. The light-gray strips represent the phase space for 
gluon splittings. If g ^qq, the quark colour factor Cp applies to part of the dipole phase 
space. The phase space reduction due to the lower energy of the produced quark is illustrated 
with the dashed strip. 



Ng{K,i). Neglecting kernel corrections and gluon splittings, this implies 



iVg(L) 



dn /^" dya{K)Ng{K), 



(24) 



giving the leading order result A^g = aA^g in Eq (p!7|a). Thus, the triangular phase space in 
Fig ^ illustrates the leading order term. 

In Fig IpD, some possible corrections up to second order are illustrated. Strips of finite width 
correspond to integrals like Eq ([2^), but with the y interval constant, resulting in terms in 
Eq (|1^) proportional to (aNg)'. A finite area, where also the k, limits are independent of L, 



represent terms proportional to {aNg 



a'N, 



At the left edge, the dark strip of width Cg/2 illustrates an effective phase space reduction due 
to kernel corrections, and the light-gray strip represents the phase-space for gluon splittings. 
To avoid drawing two strips on top of each other, the "splitting strip" is here drawn just 
inside the "Cg strip" . Together, they illustrate the first order correction —{cg-\-ni/3N^){dNgy 
in Eq (p!7|a). In a qq dipole, a corresponding strip represents —Cq^aNg)' in Eq (p!7|b). 



At the right edge, the effects of gluon sphttings is illustrated. Below the transverse scale 
of the gluon splitting, the colour factor is reduced from Nc/2 to C-p. The strips at the 
edge are changed, and due to the reduced energy of the quark, there is also a strip of 
width Cp/2 removed. These changes of strips illustrate first order corrections to the gluon 
splitting, which in itself is a first order correction, thus representing the second order term 
(ri - Cp/2)(4nf/ro3iV,)a2iVg in Eq ^ 



The colour factor reduction need not cover precisely the triangular fold below the gluon 
splitting. The border for colour suppression may be moved a finite distance, whose mean 
could be related to the average energies of the qq pair, and thus to Cp. In the generalized 
dipole equation, this mean is 3cp/4, giving the second order correction —{cpn{/2N^)a'^Ng in 
Table |1|. 

As discussed in the introduction, the top of the triangle and the surroundings of a gluon fold 
top are regions where the validity of the iterative cascade formalism is uncertain. However, 
second order corrections related to "moderately ordered gluons", emitted in these regions. 



are of some numerical importance. This is e.g. observed in [^, discussing the triangle top 
correction. 

The Cg strips of height L overshoot at the triangle top. Re-inserting a finite area, we get 
a second order correction {Cg/2)d'^Ng. Assuming the Cg strips at the edge of gluon folds to 
extend into the background triangle, we get second order corrections which precisely cancel 
the ones at the top of the gg dipole, while the top of a qq dipole is different. From the 
coefficient g2 in Table ^ we find that the dipole equation includes the (c^/2)-terms at the 
top, but not at a gluon fold, while the generalized dipole equation includes both. This is the 
main difference between the two equations, and we noted in section ^ its numerical effect on 
the anomalous dimension N' /Ng. 

Gluon splittings may also be "moderately ordered" and it is possible to have corrections 
represented by a V-shaped splitting area overshooting at the top of a gg dipole, and a 
compensating V-shaped area around the top of a gluon fold. Looking for nf dependence in 
g2, we see that these gluon splitting assumptions differ between the two dipole equations, 
but that the effects are numerically small. 

Note that the discussion here concerns possible second order corrections. One cannot con- 
clude that the equation including the most should be more appropriate. The parton cascade 
treatment of moderately ordered gluons is discussed at the end of the next subsection. 



4.2 Very Hard Gluons 

There are possible corrections related to emissions of high-fc^ gluons, not included in any 
of the evolution equations. The corrections are not numerically insignificant, but as the 
iterative cascade formalism is very uncertain for these extreme gluons, I leave the results in 
section ^ unchanged. The discussion here is mainly meant to further illustrate uncertainties 
beyond first order MLLA corrections. 



In the Dipole Cascade 

In order to take a closer look at high-fc^ corrections, we note that the emission densities off 
qq and gg dipoles are 

2(^Tr X ~\~ X X ~\~ X 

dnqq = -j^adKdy^-^-^, drigg = adndy ^ ^ ^ (25) 

In the limit kj_ <^ s, the factors dy(a;" + x^)/2 equal Kp{x)dx and i^^(x)dx, respectively. 

In the dipole evolution equations, the kernel corrections are represented by constant rapidity 
range reductions Cg and Cq. For very high-fc^ gluons, this is only approximately correct. The 
effective phase space differences, whose main contributions are at high k, have magnitude 
(index t stands for "top correction" ) 

cS'i ^ ito//'-""d.r"d!,ii±ii-i(i-c,)4 

L^oo l^Jo J~y^^^ 2 2 J 

= Cq(2-|)-ln^(2)-EA-0.23, (26) 

^ n=l ^ "- 

4^) = — + cg(2 - ^) - ln2(2) - y ^ ^ 0.40. (27) 

These top corrections modifies a2 and r2 by 

(g) (q) _ (g) Y 

Aa2 = — — ~ -0.050, Ars = — = ^ -0.042. (28) 

8 ' 4 24 ^ ^ 

Comparing with e.g. the differences in 02 between cascade formalisms, shown in Table §, we 
note that these top corrections, if trustworthy, are not numerically negligeble. 

If the only constraint on the gg dipole emission density is to reproduce the Altarelli-Parisi 
splitting function in its region of validity, the gg kernel in Eq (|25|) is not unique. We may 
generalize xl + x^ to 

xl + xl + X2[f{x2) - /(I - X2)] + ^^( Y' ^2)- (29) 

This kernel is symmetric in indices 1 and 3. The factor {k'j_/s)g does not affect Kq, provided 
g is regular. The generalized kernel in Eq (p9|) just adds a term f{x) — /(I — x) to Kq{x), 
which cancels against the corresponding term in Kq{1 — x), leaving the sum in Eq (|^) 
unaffected. (For consistency, the kernel for a mixed (qg) dipole must be modified from 

Xq + Xg to Xq + Xg + (1 — Xg)[/(x2) — /(I — X2)] + -^g{^,X2). The additional terms vanish 
when a;g — > 1, leaving Kp unchanged.) 

To examine the numerical impact of this possible kernel generelization, consider the simple 
case g = 0, /(X2) = x^2, where x is a free paramter. This modification adds x to the top 
correction cf and x/2 to the recoil coefficient 4^-* . Constraining the kernel in Eq (pOj) to be 



positive everywhere in phase space, not to overthrow the ansatz of independently emitting 
dipoles, imphes —2 < 8x < (35 + ISv^TS). With minimum %, the change in 02 is moderated 
from —0.050 to —0.034, and the change in r2 from —0.042 to —0.010, but there is no a 
priori argument for this minimum value to be appropriate. Instead letting x approach its 
maximum value ~ 10, would have drastic effects on the second order result. However, the 
possible choices in Eq (^Qf ) are constrained by data on e.g. multiplicity distributions, which 
are very well described by the cascade formalism, and four-jet characteristics in qq— i>qqgg, 
which are best described by the QCD matrix elements, but reasonably well reproduced by 
cascades |T^ . Reversing the discussion, the degrees of freedom in Eq (^) could be used to 



search for an improved simultaneous fit to multiplicity and 4-jet properties. 



In the Parton Cascade 

In the parton picture, a two-parton system is described as two independent jets, with opening 
angle 7r/2. However, in ||^ it was pointed out that this leads to an approximate treatment 
of very hard gluon emission, which has numerical consequences. To improve parton cascade 
calculations they use, for the first gluon off a qq system, the qq emission density in Eq (pSf). 



A similar well defined matrix element is not available for gg -^ ggg. If trying to include 
similar corrections in the gg case, by letting the hardest gluon be emitted in a dipole-like 
manner, rather than from two independent gluon jets, the parton formalism would inherit 
the gg dipole high-fc^ uncertainties discussed above. 



4.3 Recoil Terms 

In the dipole picture, a gg dipole takes recoils from emissions in the neighbouring dipoles, and 
thus the phase space is reduced. On average, the rapidity reduction at transverse momentum 
scale K is c^^-* /^ ~^'^ dla{l) (dashed line at left edge of Fig 0b). This gives the second order 
recoil correction —c^^^a'^Ng in Table |I[ 

A more complete set of recoil terms is illustrated in Fig ^. There, f 2 is split in two parts, 

a vr2 49 . 85 

2 6 72 ^ 72 ^ ' 

which have a direct correspondance to the quark terms ^10/^0 and wg? respectively. v\ and 
Ug come from a combination of kernel corrections with energy conservation on the emitted 
gluon, which has no correspondance in the dipole picture. 

Recoils in the parton picture, acting on an emitting gluon, corresponds in the dipole picture 
to recoils on a spectator gg dipole. Since the colour factor is N^ for gluons and A^c/2 for gg 
dipoles, this leads to the factor two difference between 4^2 and 6f^ . 

The recoil factor difference 4fio/ro — c^^ is related to the same colour factor A^c/2 as is 4vf — 
6f^ . In the emission density off a qqg system, a colour suppressed term can be represented 



10 



parton recoil terms: 



dipole recoil terms: 




dxK9.(x)ln(l-x)^4v2 



dx[K^(x) - 1 ] ln(jc) => Av^ 



N\dxK^{x)\n{l-x)^cf^ 



s(l-x) 




Cp dxE^ix) ln( 1 -x) => 4vio /''o 



dx[/^(x) - 1 ] ln(x) ^ 4v, 




2A!, 



dxi^(x) In(l-x) 
s(l-x) 



M) 




Figure 3: Comparing recoil terms, left: In the parton cascade, recoils on the emitting parton 
lead to the recoil terms 4f | and 4f lo/rg. The energy fraction x of the emitted gluon, combined 
with kernel corrections, leads to the recoil terms v\ and fg. right: In the dipole picture of 
gg — i> ggg, the squared mass of the "spectator diple" is reduced by a factor 1 — x. There 
is no corresponding spectator dipole in case of emission from a qq pair. However, a colour 
suppressed term can be represented by a dipole between the q and q, whose contribution is 
weighted with the negative factor —1/N^. This is illustrated by the dashed line in the lower 
right figure. 



by a dipole between the q and q, whose contribution is weighted with the negative factor 
— 1/N^. Recoils in the dipole picture act on this colour correction dipole. However, since it 
can not be regarded an independent emitter, full understanding of its effects requires better 
knowledge of non-perturbative QCD, e.g. questions related to so called colour reconnection. 
This is further discussed in section §. 

In the dipole model, specific assumptions about the effects of the negative dipole are made. 
Using the antenna symbol {ij) to represent a dipole between partons i and j P], the emission 
density from a qqg state can be written 



dn oc a[(qg) + (gq) - ^(qq)]- 



(31) 



The antenna {ij) has coUinear divergences in the directions of both its partons. To arrive at 
a parton cascade picture, each antenna is split in two parts, (ij) = Pij + Pji, where Pij{Pji) 
has collinear divergences in the direction of parton i{j) only, and the emission density is 
rewritten in the form 

dn ex ^[CF(Pqg + Pqg) + AT.Pgq] + a{Pgq - P,^} - — {(qq) - P^, - P^,}. (32) 

The collinear divergences cancel inside the curly brackets, and within first order MLLA 
accuracy, it is appropriate to neglect them. We note, however, that possible higher or- 



11 



der corrections from the second curly bracket, depending on the negatively weighted (qq) 
antenna, cannot be fully determined without better knowledge of non-perturbative QCD. 



5 Colour Reconnection 



In the perturbative cascade, the dipoles, or the angular ordering constraints on partons, 
associates each triplet charge with a particular anti-triplet charge. Due to the finite num- 
ber of colours, this colour topology is not unique. Interference between different topologies 
are possible, though suppressed. Furthermore, in the two most successful phenomenological 



hadronization models, a massless relativistic string [|T^ or a chain of clusters [^ connects 
triplet and anti-triplet charges. It is not necessary that this non-perturbative colour topol- 
ogy, which in the models has strong influence on the hadronic final state, agrees with the 
perturbative colour topology. 

The question of colour reconnection, where the perturbative or non-perturbative colour topol- 
ogy changes from the one defined by planar Feynman diagrams, has acquired quite some 
attention because of its possible effects on W mass measurements in e^e~ -^ W"'"W~ |^. 
To gain more insight to the problem, an observable in Z decays has been proposed P]. In 
events with a very hard gluon emission, due to which the quark and anti-quark are essentially 
collinear, colour reconnection could imply that the hard gluon forms a colour singlet system 
with some softer gluons nearby in phase space. If so, there need not be any colour flow 
connecting the qq pair and the hard gluon. This implies a multiplicity depletion at central 
rapidities along the thrust axis. 



In data from Z decays [jl5|, no effect of colour reconnection is seen. Default Monte Carlos 



are in excellent agreement with data, while a colour reconnection model [|T6| is disfavoured. 
However, in the considered events, a high-fc^ gluon is emitted, and for subsequent emissions, 
Sqq is small. According to the recoil picture where Sqq determines the phase space for colour 
factor reduction, small Sqq implies that the q and q emit coherently as an octet charge in most 
of the available phase space. In the perturbative phase, this implies that only a second gluon 
essentially collinear with the q or q could reconnect to form a singlet with the hard gluon 
in the opposite hemisphere. After such a reconnection, there is still colour flow over central 
rapidities and there is no multiplicity depletion expected. This possible qq coherence effect 



is not considered in the colour reconnection models confronted with data p, |I6|. Thus, the 
uncertainties about recoils on the negatively weighted dipole might dilute the discriminating 
power of the suggested observable. 



6 Summary 

Results from parton and dipole evolution equations, expanded one order beyond the formal 
accuracy of the modified leading log approximation, differ noticeably. This illustrates the 
theoretical uncertainty at this order. The difference between cascade pictures is larger than 
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the third order correction term calculated in the parton cascade p. Data on Ng/N^ |^, || 
seem to favour the smaller result in the dipole cascade. 

One important source of uncertainties are "moderately ordered gluons", i.e. two gluons with 
similar k±, or a first hard gluon with k± close to the kinematical limit. In particular, 
the emission density for very hard gluons off gluonic systems is not fully determined by 
the constraint on soft gluon emission given by the Altarelli-Parisi splitting function. This 
gives an opportunity to search for a emission kernel which improves the simultaneous fit to 
multiplicity properties (very well described by cascades) and four-jet characteristics (better 
described by matrix elements). 

Part of the recoil effects are truly non-perturbative, related to negative colour suppressed 
terms which can not be represented by independently emitting sources. There are possible 
relations between these recoil effects and the question of colour reconnection. 

In this paper a generalized dipole evolution equation is presented, which includes energy 
conservation effects in a more general way than has been done before in the dipole formal- 
ism 10. This equation is suitable to calculate higher order corrections to all multiplicity 
moments, not only average multiplicity, in the dipole cascade picture. 
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A The Generalized Dipole Evolution Equation 

In this appendix, I derive a generalized dipole evolution equation alternative to the one 
presented in ref [|10|. Since the dipole cascade is ordered in k_j_, I introduce a generating 



function P(L, k), where L and k are logarithms of the dipole energy and maximum allowed 
k±, respectively. 



Gluon Dynamics 

In gluon dynamics, based on dngg in Eq (pSD , the generating function for dipoles satisfies 

rl{L,K) = / dK'2a{K')S{L,K,K') dy^^Pg(Li2,«:')^g(L23,«:')^g(i^i3,«:') 

+S{L,K,K,)Vl{L,K,), (33) 

where ymax is defined in Eq (^3]), and the Sudakov form factor 



S{L, K, k') = exp 



d."2a(.") r^ dy^ 



"^niax 



(34) 
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ensures that the emitted gluon is the one with highest k_j_. An infrared cut-off kj_cut is 
introduced, with Kc = lii(^icut)- Multiplying by S~^{L, k, Kc) and taking the derivative 
w.r.t. K leads to, with Vg = Vg{L, k), Vg-- = Vg{Lij, n), 



&P /-J/max 






-p -p -p — "P ^ 

' gl2' 623' S13 ' S 



When K, <^ L, this has a solution linear in L. Parametrizing the n dependence as 

ln(Pg(L, k)) = Agin) + {L- K- Cg)B{K) 

gives a pair of coupled evolution equations 

A'^-B = aexp{Ag-CgB)J2Rlf^B^, 



n=l 



B' 



a 



[exp{Ag-CgB) -1], 



where 



1 /"l r\T 1 1 

C being the Riemann zeta-function. 



(35) 

(36) 

(37) 
(38) 

(39) 



qq Systems 

Gluon emission off a qq system creates a qg and a gq dipole. In addition to these, a colour- 
correction term can be represented by a dipole between q and q, whose emission density is 
weighted with —1/N^. Following the arguments in |jT^, its only effect can be assumed to be 
a reduction of the colour factor in "mixed" qg dipoles from Nc/2 to Cp, in a rapidity range 
Yq^ near the quark. For simplicity, quark masses are here neglected. 

With Vm the generating function for the multiplicity distribution in a mixed dipole, we get 
evolution equations in the form 



dn 
dVmJY) 



OC^-r^ rymax T -I- T 

^a{K) / dy^l^ [Vrm2{Yq)Vm2AY^} - PJ , 



(40) 



a[K] 



Vn 



dy 



xf + xl 



ymax 

1 

c 



V 

'PmuiX )^g23^^~~ 



VmiY) 



SR 



X 



e(r-y„ 



y) 



(41) 



Here Y' represents the rapidity range of colour suppression in the quark direction, after the 
gluon is emitted from a mixed dipole. The original gluon of the mixed dipole is assumed to 
be attached to a gg dipole with logarithmic scale Lr, which due to recoils is reduced to Lri. 
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This equation contains a set of more or less free parameters, but reasonable assumptions 
implies 



ln(Pq(L,/€)) = Ac^{k) 



2Cf 

1^ 



{L — K — Cq)i3(/€) 



(42) 



ln{VmiL,n,Y)) = ^[\niV^(L,K)) + lniV,{L,K)) + ^{L-2Y)BiK)]. (43) 



This adds to the gluon dynamics result, Eq (^) and (|38|), a quark evolution equation 

TV 
"-Aq' = B + a exp(^g - c^B) J2 R^n^^". 

' n=l 



2Ci 



(44) 



where R^' depend on the assumtions made on Yq and 1^ in Eq (1^) . In [|T^ , two reasonable 
possibilities are examined. One is Fq + 1^ = ln(s) = L, which implies that all i?^^^ are 0. 
The other is Yq + Ig- = In(sqq) = L13, which implies 

1 _i /-i Hr 1 1 

^^'^ = ^Sm^^ i -[' + (' - ")'] ^""^^^ - "^ = (]^)"[2C(n + 1) - 1 - ^]. (45) 



There are some indications from experimental data in favour of this second choice ^ and 
we will limit ourselves to this in the present paper. 



Pair Production 



To include pair production, g— >qq, (still neglecting quark masses) we add a gluon splitting 
term 



rifayK,) 



f dX2(x2 + xl)[Vml2{Y^)Vml3{YTd - Pg 

J2k I /\/s 



2Nc J2kj_/y^ 



(46) 



to the gluon evolution equation, Eq (pSj), and a similar term to the evolution equation for a 
mixed dipole, Eq (0). With reasonable assumptions about recoil effects when the gluon in 
a mixed dipole sphts into a qq pair, the only modification induced by gluon splittings is an 
additional term 



-a 



3iV, 



, . . r2C/p 

exp(>lq — Ag — [— — Cq — C| 



.}B) y: b-r^^^ - 1 



n=0 



(47) 



on the right-hand side of Eq (^Tf ). Except for Rq = 1, R^^^ depend on the assumed 
magnitude of the colour reduction phase space Yq + Y^ near the new qq pair. The intuitive 
assumption Yq + Y^ = In(sqq) = L23 implies 



1 L^ 



2n\Jo 



/ dx{x'^ + (1 — x) 
Jo 



;i + — )ln(x) + ln(l-a;) 



(48) 
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Combined Result 

Neglecting possible corrections related to emissions with k± close to the kinematical limit, 
discussed in section ^^ , the solutions to ln(P(L, k)), Eq ( ^6l) and (H), are valid also when 
K, is close to L. This implies Ag{L — Cg) = ln(Pg(L)) and Aq{L — Cq) = ln('Pq(L)), where 
Vg{L) and Vq{L) represent unbiased dipoles, i.e. dipoles with only kinematical constraints 
on k±. Thus, Eqs (^, (|S|), (|§ and (g^) combine to Eqs (|IB|). 
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